Could Z6(10610) be a B*B molecular state? 
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Assuming tlie newly observed structure ^6(10610) as a bottomonium-like molecular state B*B, 
we calculate its mass in the framework of QCD sum rules. The numerical result is 10.54 ±0.22 GeV 
for B* B, which coincide with the mass of 2'i,(10610). This consolidates the statement made by Belle 
Collaboration that the ^(,(10610) resonance could be a B* B molecular state. 
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^ ; I. INTRODUCTION 

^ ■ Very recently, Belle Collaboration observed two narrow structures ^6(10610) and .^b(10650) in the 
7r='=T(n5') (n — 1,2,3) and n'^hbimP) {m = 1,2) mass spectra that are produced in association with a 
■ single charged pion in T{5S) decays [l|. The measured masses of the two structures are 10608.4 ± 2.0 MeV 
and 10653.2 ± 1.5 MeV, respectively. Experimental analysis favors quantum numbers of J = 1 for both 
states. As Zf,(10610) and Zf,(10650) are charged, they can not be simple bb combinations. The measured 
rS i' masses of these two new states exceed by only a few MeV the thresholds for the open beauty channels 
ji^ B*B and B*B* . They could be interpreted as molecular states and determined by the strong interaction 
dynamics of B*B and B* B* meson pairs 

The concepts of molecular states were put forward long ago in f3|. Some of the exotic X, Y, and Z 
resonances have been described as possible charmonium-like molecular candidates in the literatures since 
their masses are very close to the meson- meson thresholds. Explicitly, it is interpreted Z+(4430) as a 
D*Di molecular state 0^(3930) as a D*D* molecular state jJH], y(4140) as a DJZ)* molecular state 
, y(4260) as a XcP° llO| or an ujXci molecular state [11|, X(4350) as a D*D*q molecular state |12l.ll3j. 
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y(4274) as a I?s£'so(2317) molecular state [lj| etc.. If molecular states can be confirmed, QCD will be 
further testified and then one will understand QCD low-energy behaviors more deeply. 

The newly observed Zh resonances may open a new window to study molecular states in the 
bottomonium-like family. Therefore, it is interesting to investigate whether they could be bottomonium- 
like molecular candidates. The quantitative description of their properties like masses are helpful for 
understanding their structures. Unfortunately, quarks are confined inside hadrons in the real world, and 
the strong interaction dynamics of these states are governed by nonperturbative QCD effect completely. 
^ In this work, by assuming Zb(10610) as a B*B molecular state, we calculate the mass of this resonance 



in the framework of QCD sum rule (QCDSR) method [15|, which is a nonperturbative formulation fi rmly 
rooted in QCD basic theory and has been used to study some charmonium-like molecular states 16l- 18|. 



It is not so straightforward from the meson-meson configuration of fields to construct a B*B* current with 
a quantum number of = 1+. The Zf,(10650) will not be discussed here. Our final numerical result 
10.54 ± 0.22 GeV for B*B agrees well with the experimental data of Z'b(10610), while the masses of the 

which are 
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= 1+ bbqq tetraquark states were found to be around 10.1 ~ 10.3 GeV in QCDSR 
lower than the measured values of Zi, states. The present work thus favors that the Zb(10610) resonance 
could be a B*B molecular state rather than a bbqq tetraquark state. 

The rest of the paper is organized as three parts. We discuss QCD sum rules for the molecular state 
in Sec. |lTl where the phenomenological representation and the operator product expansion (OPE) contri- 
bution up to dimension six operators for the two-point correlator are derived. The numerical analysis is 
made in Sec. IIIII The mass of the B*B molecular state is extracted out and found to coincide with the 



experimental value of Zb(10610) resonance. The Sec. IIVI is a short summary and outlook. 
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II. MOLECULAR STATE QCD SUM RULES 

The starting point of the QCD sum rule method is to construct the interpolating current properly 



and then write down the correlator (for reviews see (2l|-|2J| and references therein). The molecular state 



currents are built up with the color-singlet currents for their composed hadrons. As for B* B, the current 
is constructed as 

fs-B ^ iQel''be)(bfi-/5qf), (1) 

where q indicates the light quark and the subscript e and / are color indices. Note that the current is local 
and the four field operators act at the same space-time point. It is a limitation inherent in the QCDSR 
disposal of the hadrons since the bound states are not point particles in a rigorous manner. The current 
is different from that of tetraquark state which is diquark-antidiquark configuration of fields. These two 
types of currents can be related to each other by Fiertz rearrangement and differ by color and Dirac factors 



25| . It will have a maximum overlap for the molecular state with the meson- meson type current. The 



masses of the bbqq tetraquark states were calculated in QCDSR and found to be around 10.1 ^ 10.3 GeV 



in [19j, and 10144 ± 106 MeV in [20[. Both of the predicted values for tetraquark states are lower than the 



Zh mass. Therefore, it's hard to accommodate Zf, as tetraquark states. 
The two-point correlator is defined as 



Lorentz covariance implies that it can be generally parameterized as 

nA'-(g2) = _ gA.-)n(i)(g2) + ^n(o)(g2). (3) 

The term proportional to g^y will be chosen to extract the mass sum rule. In phenomenology, Il^^\q^) 
can be expressed as 

nu.,,^,.-^ + ir.,£Hl5^, (4) 



B*B 



where Mg.g denotes the mass of the B*B resonance, sq is the threshold parameter, and A'^-'^^ gives the 
coupling of the current to the hadron q\B* B) = A^^^e''. In the OPE side, IV'^^\q^) can be written as 



nW(g^)=/ dsP ^, (5) 

J Ami * ~ ^ 

where the spectral density is p^^^{s) = ilmll^"'^) (s). Applying quark-hadron duality and making a Borel 
transform, we have the sum rule from Eqs. @ and ([5]) 
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with AP the Borel parameter. Making the derivative in terms of M"^ to the sum rule and then dividing 
by itself, we have the mass of the B* B state 

J Ami J Ami 

For the OPE calculations, we work at leading order in Ug and considers condensates up to dimension six. 



with the similar techniques developed in [26|. To keep the heavy-quark mass finite, the momentum-space 
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expression for the heavy-quark propagator and the expressions with two and three gluons attached are 
used fi?!. The hght-quark part of the correlation function is calculated in the coordinate space and then 
Fourier-transformed to the momentum space in D dimension. The resulting light-quark part is combined 
with the heavy-quark part before it is dimensionally regularized at 13 = 4. The spectral density can be 
written as 

pOPE(^) ^ pPert(5) + p(99>(,) +p{??>^(,) +pto-G,>(,) +^{9=G^>(,) +^(S^G^>(,)^ (8) 

where pP^^^, p^^'^\ p^^'^'^ , p'^9q<^-Gq) ^ p{g G ) ^ ^^^^ p(g G^) ^^^.g -j-j^g perturbative, quark condensate, four- 
quark condensate, mixed condensate, two-gluon condensate, and three-gluon condensate spectral densities, 
respectively. They are 



Pmir. 



with r(mh,s) — (a + li)m^ — afis. The integration limits are given by amin ~ ~ -^/l — Am'^/ /2, 

ctmax = ^H-a/I — 4m^/s^ /2 and Pmm — arn^/ {sa—ml). We have applied {qqqq) = K{qq)'^ to estimate the 
four-quark condensate, where the parameter k is introduced to account for deviation from the factorization 
hypothesis 22 1. 



III. NUMERICAL ANALYSIS 



In this section, the Eq.Q will be numerically analyzed. The b quark mass is taken as rrib = 4.20j^Q j!)y GeV 



-0.07 

28|. The condensates are (qq) = -(0.23 ± 0.03)3 QgV^ {gqa ■ Gq) = ml (qq), = 0.8 GeV^ (g^G^) 



0.88 GeV'', and {g^G^) = 0.045 GeV'^ [221, |26|. In the QCDSR approach, there is approximation in the 
OPE of the correlation function, and there is a very complicated and largely unknown structure of the 
hadronic dispersion integral in the phenomenological side. Therefore, the match of the two sides is not 
independent of M^. One expects that there exists a range of M^, in which the two sides have a good 
overlap and the sum rule can work well. In practice, one can analyse the convergence in the OPE side and 
the pole contribution dominance in the phenomenological side to determine the allowed Borel window: on 
one hand, the lower constraint for is obtained by the consideration that the perturbative contributions 
should be larger than the condensate contributions, so that the convergence of the OPE is under control 
and the higher dimension terms can be safely ignored; on the other hand, the upper limit for is obtained 
by the restriction that the pole contributions should be larger than the continuum state contributions, so 
as to guarantee that the contributions from high resonance states and continuum states remains a small 
part in the phenomenological side. Meanwhile, the threshold parameter is not completely arbitrary 
but characterizes the beginning of the continuum states. The energy gap between the groundstate and 
the first excitation is around 500 MeV in many cases of nucleons or charmonium-likc states. Whereas, 
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FIG. 1: The phenomenological contribution in sum rule ([6)1 for ^/so — 11.4 GeV. The sohd hue is the relative pole 
contribution (the pole contribution divided by the total, pole plus continuum contribution) as a function of 
and the dashed line is the relative continuum contribution. 

this may not be always straightforwardly generalized to the bottomonium-like states since their absolute 
masses are much larger than masses of nucleons or charmonium-like states. One should consider whether 
there is a proper region of Borel parameter Af ^ for a fixing ^/sq, or whether one could find a working 
window for the sum rule. On all accounts, it is expected that the two sides have a good overlap in the 
determined work window and information on the resonance can be safely extracted. 

Firstly, we take the factorization hypothesis of the four-quark condensate {qqqq) — K{qq)'^ with k — 1. 
The comparison between pole and continuum contributions of sum rule © as a function of the Borel 
parameter for the threshold value ^/so = 11.4 GeV is shown in FIG. 1. Its OPE convergence by 
comparing the perturbative, quark condensate, four-quark condensate, mixed condensate, two-gluon con- 
densate, and three-gluon condensate contributions as a function of is shown in FIG. 2. The ratio 
of perturbative contributions to the total OPE contributions at = 9.0 GeV^ is nearly 58% and it 
increases with M^. Thus the perturbative contributions will dominate in the total OPE contributions 
when > 9.0 GeV^. On the other hand, the relative pole contribution is approximate to 50% at 
— 10.5 GeV^ and it decreases with AP. In order to guarantee that the pole contribution can dominate 
in the total contributions, we have the value AP < 10.5 GeV^. Thus, the range of A/P for B* B is taken 
as = 9.0 ~ 10.5 GeV^ for = 11.4 GeV. Similarly, the proper range of A^P is 9.0 - 9.8 GeV^ for 
y/so = 11.2 GeV, and 9.0 ^ 11.2 GeV^ for y/so — 11.6 GeV. We see that the corresponding Borel parameter 
range is A/P = 9.0 - 9.1 GeV^ for ^ = 11.0 GeV, and AP = 9.0 - 9.4 GeV^ for ^ = 11.1 GeV, which 
are very narrow as working windows. Therefore, the threshold parameter y'so is taken as 11.2 ^ 11.6 GeV. 
The mass of B*B is numerically calculated to be 10.56 ± 0.18 GeV and shown in FIG. 3. 

The coupling constant A*-^^ between current and the particle is calculated from Eq.® in the same 
working windows. We arrive at A'^-* — 0.27 ± 0.07 GeV^. The equivalent quantity of the tetraquark 
current is 0.09 ^ 0.11 GcV^ in [23|- The coupling constant A*^^-' is roughly three times as large as the one 
of the tetraquark current and the meson-meson molecular current has a larger overlap with the Z\, state 
in comparison to the diquark-antidiquark tetraquark current. 

To investigate the effect of the factorization breaking, we assume that (qqqq) ~ i^{qq)'^ with k — 2. 
From the similar analysis process, the corresponding working windows are taken as: A/P — 9.0 ^ 9.9 GeV^ 
for ^ = 11.2 GeV, AP = 9.0 - 10.6 GeV^ for ^/I^ = 11.4 GeV, and = 9.0 - 11.4 GeV^ for 
^/Jq = 11.6 GeV. We extract the mass value 10.52 ± 0.21 GeV. Finally, we average two results for 
K = 1, 2 and arrive at the mass value 10.54 ± 0.22 GeV for B*B, which agrees with the experimental value 
10608.4 ± 2.0 MeV for Zb(10610). 
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FIG. 2: The OPE contribution in sum rule ((6]) for ^/sQ = 11.4 GeV. The OPE convergence is shown by comparing 
the perturbative, quark condensate, four-quark condensate, mixed condensate, two-gluon condensate, and three- 
gluon condensate contributions. 




FIG. 3: The mass of the B* B molecular state as a function of from sum rule ((TJ. The continuum thresholds 
are taken as ^ = 11.2 ~ 11.6 GeV. The ranges of is 9.0 ~ 9.8 GeV^ for ^ = 11.2 GeV, 9.0 ~ 10.5 GeV^ 
for = 11.4 GeV, and 9.0 ~ 11.2 GeV^ for ^ = 11.6 GeV. 



IV. SUMMARY AND OUTLOOK 



By assuming Zb(10610) as a B*B molecular state, the QCD sum rule method has been applied to 
calculate the mass of the resonance. Our numerical result is 10.54 ± 0.22 GeV for B*B. It is compatible 
with the newly measured experimental data of ^(,(10610) by Belle Collaboration, which supports the 
statement that Zb(10610) resonance could be a B*B molecular state. It is expected that this work is 
helpful for understanding the structure of Zb(10610). For the newly observed structure ^(,(10650), one 
could consider how to construct a B* B* molecular state current with a quantum number of — 1^ from 
the meson-meson configuration of fields. For further work, one needs to take into account other dynamical 
analysis to identify the structures of Zt hadrons. 
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